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We study the attractive Hubbard model with mass imbalance to clarify low temperature properties of the 
fermionic mixtures in the optical lattice. By combining dynamical mean-field theory with the continuous- 
time quantum Monte Carlo simulation, we discuss the competition between the superfluid and density wave 
states at half filling. By calculating the energy and the order parameter for each state, we clarify that the 
coexisting (supersolid) state, where the density wave and superfiuid states are degenerate, is realized in the 
system. We then determine the phase diagram at finite temperatures. 
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Superfluid state in ultracold fermionic systems has 
attracted considerable interest since the successful re- 
alization of the Bose-Einstein condensation in ^Li2 
molecules.^' Due to the high controllability in the sys- 
tem, remarkable phenomena have been observed such as 
the BCS-BEC crossover^-'*) and the superfluid state in 
the spin-imbalanced system,^' where Cooper pairs are 
composed of ions with distinct hyperfine states. Recently, 
the fermionic mixtures with distinct ions, eg. ^Li and 
^''K, have experimentally been realized,^' which stimu- 
lates further theoretical investigations on the superfluid 
states in the mass imbalanced system. ^"^''^ 

One of interesting questions in such a mass imbalanced 
system is how the superfluid state is realized when the 
lattice potential is loaded, so-called, an optical lattice. 
In the lattice system, the density wave (DW) state is 
naively expected, in addition to the SF state, since less 
mobile fermions tend to crystallize in the lattice, particu- 
larly, at half filling. It is desired to systematically discuss 
how the SF state competes or coexists with the DW state 
in the optical lattice system. This topic is closely related 
to an important issue in condensed matter physics, so- 
called, the supersolid state, since the DW state can 
be regarded as a sort of the solid state. Therefore, the 
optical lattice system with the mass imbalance should be 
providing an ideal stage for the studies of the supersolid 
state in fermionic systems. 

Motivated by this, we study low temperature proper- 
ties in the fermionic mixture in the optical lattice, com- 
bining dynamical mean-field theory(DMFT)^^^ with the 
continuous-time quantum Monte Carlo (CTQMC) simu- 
lations.^®' By calculating the order parameters for the 
DW and SF states, we determine the phase diagram at 
finite temperatures and clarify how the coexisting state 
is stabilized against the mass imbalance. 

In this paper, we consider the following attractive Hub- 
bard model with different masses, as 

H = Uictcj^ + h.c.) - U^n,^n,i (1) 
where denotes nearest neighbor site, cl^{cia) is the 



creation (annihilation) operator of a fermion at the ith 
site with spin a{=1',X) and Uia- = cl^Cia. U{> 0) is the 
attractive interaction and ta is the hopping amplitude 
for the fermion with spin a, where the effect of the mass 
imbalance is taken into account. 

We examine low temperature properties of this model 
by means of DMFT^^^ which maps the lattice model to 
the problem of a single-impurity connected dynamically 
to a "heat bath". The Green's function is obtained via 
the self-consistency condition imposed on this impurity 
problem. When the DW and SF instabilities are equally 
treated in the DMFT framework, the self-consistency 
equation for the sublattice a[= A, B] is given as^^-* 



Goa{z] 



zero + l^cr, 



-ifGa(z)f, 



(2) 



where ctq is the identity matrix, dz is the z-component 
of the Pauli matrix. ^ is the chemical potential and 
T = diag (I?!, — Uj,), where is the half bandwidth 
for the bare band with spin a. Goa{z) and Gq(z) are the 
noninteracting Green function for the effective impurity 
model and the local Green function for the sublattice a, 
which are represented in the Nambu formalism. 

There are various methods to solve the effective im- 
purity problem. To study how the competition between 
the DW and SF states in the mass imbalanced system, an 
unbiased and accurate numerical solver is necessary, such 
as the exact diagonalization or the numerical renormal- 
ization group. A particularly powerful method for explor- 
ing finite temperature properties is the hybridization- 
expansion CTQMC method,^®' which enables us to 
study the attractive Hubbard model both in the weak- 
and strong-coupling regimes. ^^-^ In the paper, by vary- 
ing the ratio of the bandwidths r = D^/ with a fixed 
= 1 (energy unit), we proceed to discuss how the 
mass imbalance affects low temperature properties. 

In the mass balanced system (r = 1), low-energy prop- 
erties have been studied in one dimension, ^^"^^^ two di- 
mensions'^^' '^^^ and infinite dimensions. It is known 
that the DW and SF states are degenerate on the bipar- 
tite lattice in two and higher dimensions at half filling. 
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Each order parameter can be defined by 



A 



DW 



(3) 



(4) 



In this case, the supersolid state with both order param- 
eters is possible to be reahzed. On the other hand, in the 
imbalanced hmit (r = 0), the system is reduced to the 
spinless Fahcov-Kimbah model with mobile and localized 
fermions,'*^"^^-' where the ground state is the genuine DW 
state. Therefore, it is necessary to clarify how the intro- 
duction of the mass imbalance lifts the degeneracy of 
these two states. 

To clarify this, we first discuss the stability of each 
state in the system. We calculate the order parame- 
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Fig. 1. Order parameters for density wave (a) and superfluid (b) 
states as a function of the interaction strength at T = 0.05. 



ters Aow and Asf under the conditions Asf — and 
Adw — 0, respectively. The obtained results for the DW 
and SF states at T = 0.05 are shown in Fig. 1. In the 
mass balanced case (r = 1), the system has the SU{2) 
symmetry and thereby we find that these two order pa- 
rameters are identified within the numerical accuracy. 
In the case, the nature of the phase transition is well- 
known.'^'^-' In the noninteracting case [/ = 0, a normal 
metallic state is realized. Increasing the interaction be- 
yond Uci{— 0.95) the order parameters are induced, and 
the phase transition occurs to the DW or SF state. Fur- 



ther increase in the interaction drives the system to the 
normal metallic state at Uc2{= 5.0). 

The introduction of the mass imbalance leads to dif- 
ferent behavior where A^w > Asf, as shown in Fig. 1. 
When r = 0.4, the critical interactions for both states 
are deduced as Uci{DW) = 0.83, U^2{DW) = 2.6, 
UaiiSF) = 0.96, and Uc2iSF) = 2.0, by examining the 
critical behavior A ~ (|C/ — Uc\/Uc)^ with the exponent 
(3 = 1/2. These results mean that the DW state is stable 
in the larger parameter space while the SF state becomes 
unstable. In fact, when r < 0.3, the SF solution no longer 
exists and the genuine DW state is realized in the inter- 
mediate coupling region. 

By performing similar calculations, we obtain the 
phase diagram at T = 0.05, as shown in Fig. 2. In the 
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Fig. 2. Phase diagram of the half-filled attractive Hubbard 
model at T = 0.05, where the lines are guides to eyes. Solid lines 
represent the phase boundaries and a dashed line represents the 
boundary where the unstable SF solution disappears. Open (solid) 
squares indicate the coexisting (genuine DW) state. 



weak and strong coupling regions, the normal metallic 
state is realized due to thermal fluctuations. When r — \, 
the coexisting state appears in the region (J7ci < U < 
/7c2), as discussed above. Introducing the mass imbal- 
ance on the system with J7 = 2, the order parameters for 
the DW and SF states are monotonically decreased and 
at last disappear at the critical points r^^('^ 0.15) and 

^SF 



(~ 0.42), as shown in Fig. 3 (a). To study how the 
degeneracy of these two states is lifted, we also calculate 
lergy for eacl 
£:,7,i9) with 



the internal energy for each state i^DW — + Eu and 



SF 



= Ef 



E 



DW 
K 



E 



2 



/3 



dTGAa{T)GBa{-T), (5) 
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dT 



E 



Da 

2 
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where G(r) and F{t) are the normal and anomalous 
Green's functions. The results are shown in Fig. 3 (b). 
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Fig. 4. Order parameters for the density wave (a) and superfiuid 
(b) states as a function of the temperature T in the system with r = 
0.2,0.4,0.6,0.8 and 1.0, which are obtained under the conditions 
^SF = and ^DW = 0, respectively. 



Fig. 3. Order parameters (a) and internal energies (b) as a func- 
tion of the ratio r in the optical lattice with U = 2 at T = 0.05. 



It is found that when Tc < r < 1^ the internal ener- 
gies are identified within the numerical accuracy, where 
Tc ~ 0.8 ~ 0.9. Therefore, we can say that the coex- 
isting state is realized in the region. Below r — r^^ the 
degeneracy of two states is lifted, where the genuine DW 
state is realized and the SF state becomes unstable. Fur- 
ther decrease in the ratio r yields the second-order phase 
transition to the normal state at the critical value , 
where the order parameter vanishes and the curve of the 
energy for the DW state merges with the paramagnetic 
one. In the limit r = 0, the normal metallic state appears 
at T = 0.05, which is consistent with the result for the 
Falicov-Kimball model. 

Fig. 4 shows the temperature dependence of the or- 
der parameters for the SF and DW states in the system 
with U ~ 2, which are obtained under the conditions 
^DW = and Asf = 0, respectively. We find that when 

< r < 1, the magnitudes of these quantities are iden- 
tified. This implies that the coexisting state is realized 
below the critical temperature. On the other hand, in the 
small r case, the critical temperatures for the DW and 
SF states are difi^erent from each other, e.g. T^^ — 0.065 



and T^^ = 0.076 when r = 0.6. This implies that the 
SF state is unstable, and the genuine DW state is realized 
below T^^. 

We wish to discuss ground state properties in the 
strong coupling limit. To this end, we first use the 
particle-hole transformation'^"' as ai-^ = Cif and a^^ — 
(— l)*c|^j^, and map the negative-[/ Hubbard model to the 
positive-?/ Hubbard model. Then, the effective Hamilto- 
nian at half filling is given by the XXZ quantum spin 
s = 1/2 model as, 



H 



(8) 



where Si = 5 ^L'^a/J^*/?' ^ = (*t " *4.)^/2it^i and 
J = At-ft^/U. In the mass balanced case (r = 1), the 
effective model is reduced to the isotropic Heisenberg 
model. The ground state is the antiferromagnetically or- 
dered state, where the direction of the ordered moment is 
arbitrary. Namely, when the ordered moments are along 
the z axis (in the x — y plane) , the DW (SF) state is re- 
alized in the original attractive Hubbard model. On the 
other hand, the mass imbalance yields the anisotropy 6 
in the spin Hamiltonian, where the ordered moments are 
fixed along the z axis. This implies that in the strong cou- 
pling limit, the coexisting state is realized only at r = 1 
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and the genuine DW state is realized in general. There- 
fore, we can say that in the attractive Hubbard model 
with different masses, the phase boundary between the 
DW and coexisting states approaches the r = 1 axis when 
U increases. The coexisting state is then realized only in 
the intermediate coupling region. It is expected that the 
hole doping induces the genuine SF state (phase separa- 
tion) in the coexisting (genuine DW) region of the phase 
diagram at half filling, which is consistent with the re- 
sults obtained from DMFT with the exact diagonaliza- 
tion.i9) 

Before closing the paper, we would like to comment 

on the realization of the supersolid state in the fermionic 
mixtures on the optical lattice. We have confirmed that 
the coexisting state appears in the phase diagram. It is 
expected that such interesting behavior appears in the 
three dimensional optical lattice with the mass imbal- 
ance. However, one may consider that it is difficult to re- 
alize the supersolid state since SF fluctuations are hard 
to be controlled in the system. Nevertheless, there are 
some possibilities to realize the supersolid state in the 
fermionic system. One of them is the introduction of 
frustration to the DW state realized at low temperatures 
since it tends to destabilize the DW state and to enhance 
SF fluctuations. Another possibility is the introduction 
of the lattice potential to the superfluid state realized in 
the fermionic mixtures. By tuning DW and SF fluctu- 
ations by means of corresponding parameters properly, 
the supersolid state should be realized in the fermionic 
mixtures on the optical lattice experimentally. 

In summary, we have investigated low temperature 
properties of the fermionic mixtures in the optical lat- 
tice, which should be described by the attractive Hub- 
bard model with different masses. By combining DMFT 
with the strong-coupling version of the CTQMC method, 
we have studied half-filling properties carefully to clarify 
that the coexisting (supersolid) state, where the DW and 
SF states are degenerate, is realized. By performing sys- 
tematic calculations, we have obtained a rich phase dia- 
gram in the mass imbalanced attractive Hubbard model. 
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